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Abstract — We derive a closed-form expression for the or- 
thogonal polynomials associated with the general lognormal 
density. The result can be utilized to construct easily computable 
approximations for probability density function of a product 
of random variables, when the considered variates are either 
independent or correlated. As an example, we have calculated 
the approximative distribution for the product of Nakagami-m 
variables. Simulations indicate that accuracy of the proposed 
approximation is good with small cross-correlations under light 
fading condition. 

Index Terms — Product of random variables, central limit theo- 
rem, lognormal distribution, orthogonal polynomials, Nakagami- 



I. Introduction 

STATISTICAL properties of products of Random Variables 
(RVs) are essential in performance analysis of contempo- 
rary wireless communication systems. For example, the fading 
amplitude of multi-hop relaying systems follows the distribu- 
tion of the product of Nakagami-m RVs |1|. In addition, the 
cascaded-keyhole channel can be modeled using the product 
of individual keyhole channels Q . 

The exact Probability Density Function (PDF) and Cumula- 
tive Distribution Function (CDF) of the product of independent 
Beta, Gamma and Gaussian RVs can be represented in terms 
of the Meijer-G function (3\. A more general framework 
involving the Fox H-function was proposed in [4J for the 
distribution of product of almost any non-negative independent 
RVs. Although being theoretically interesting, these G- or H- 
function representations are difficult to evaluate on a low- 
complexity computing platform, where the calculation of the 
product statistics is required [5]. The numerical calculations 
of G- and H-functions require a numerical solution of contour 
integrals or a construction of look-up tables that cover all 
parameter combinations. Both approaches set stringent re- 
quirements for the considered platform either in the computing 
capability or the memory capacity. Recently, Ahmed et al. ||6l 
employed the Mellin transform and the residue theorem to 
derive the PDF of a product of i.i.d. Nakagami-?7i RVs as 
an infinite series. By exploiting the structure of the Mellin 
transforms, the authors in |7|, |8| proposed approximations 
to the PDFs and CDFs of products of independent Rayleigh, 
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Gamma, Nakagami-m and Gaussian RVs. However, when the 
RVs in product are correlated, exact distributions are known 
only in some special cases (e.g. bivariate Nakagami-m in [9|). 
For arbitrary number of RVs with generic correlations, none 
of the aforementioned methodologies are suitable to give a 
tractable solution. This is the main issue to be addressed in 
this paper 

Motivated by the Central Limit Theorem (CLT), the approx- 
imation for the product of RVs can be constructed by using a 
lognormal density and its associated orthogonal polynomials. 
We apply the approach of |10| where the first few moments 
of the product and the approximated density function are set 
equal. The resulting approximative distribution involves only 
a finite sum of polynomials and the elementary lognormal 
density function. Therefore, implementing the proposed ex- 
pressions on practical computing platforms is straightforward 
compared with the G- or H-function representations. We note 
that the proposed framework is suitable for cases of both 
independent and correlated variables provided that moments 
of a product can be computed in closed-form. 

In this paper we deduce a closed-form expression for the 
orthogonal polynomials associated with a general lognormal 
density. The result is obtained by taking advantage of the 
determinant representation of orthogonal polynomials. To give 
a specific example, we derive an approximative PDF and 
CDF for the product of both independent and correlated 
Nakagami-771 RVs. Numerical results are compared with sim- 
ulations in terms of the Complementary CDF (CCDF). The 
overall approximation accuracy is measured by the Mean 
Square Error (MSE) between the approximative CDF and the 
empirical CDF. The orthogonal polynomials associated with 
a standard lognormal densitjo were first developed in [11], 
which is a special case of our result. It is worth of noticing 
that our results are not straightforward extensions of the results 
presented in 1 11 1 but require new methodological elements. 

II. Orthogonal Polynomials Associated with 
General Lognormal Density 

Assume that y is a Gaussian RV with mean /i and variance 
(7^. Then X ~ e^ follows the lognormal distribution 

/LN(a;) = -^==e ^"^ , a;e(0,oo) (1) 

'When the corresponding Gaussian distribution is zero mean and unit 
variance. 



and the ^-th moment, Ui (i € N), of X is given by 
Let 7r„(x) be the n-th degree polynomial 






(2) 



(3) 



where c„_„ j^ 0. The polynomials {7r„(x)} are said to be 
orthogonal with respect to the general lognormal density if 



TTj{x)TTk{x)fi^j^{x)dx = hj6jk, j, fc e N, (4) 



where hj — J2i=o ^k=o '^j.i^^j.k'^i+k is a normalization factor 
The symbol Sjk is the Kronecker delta symbol, which is 
defined as Sjk = 1 if j ~ k and zero otherwise. 

Due to (|2|i any arbitrary moment of fLN{x) exists. Thus, 
functions x' (i G N) belong to the space of square inte- 
grable functions with respect to the weight function fLN{x)- 
There exists a unique set of orthogonal polynomials {7r„(a;)} 
which admits the explicit determinant representation 7r„(a;) = 
(A„)"iA„(a;) lUll Th. 2.1.1], where A„(a;) denotes the n-th 
degree polynomial 



A„(x) = 



^0 



Vn-l 
1 



i^n-1 






(5) 



and the constant A„ is formed by deleting the last row and 

column from A„(x), i.e. A„ = |i^i+ili.j^o.-- ,«-i ^i* ^o "" 
1. The determinant (|5]l can be expanded using the cofactors 
with respect to the last row. Then 7r„(a;) becomes 



" A 



(6) 



fe=0 



where the cof actor A„ ^ with respect to x^ is obtained by 
removing the (fc + l)-th column and the last row from (|5]l, 
Ao.o = and Aio = 1. After comparing (|6]l with (|3}, we 

obtain 

, A., I, 

(7) 






Let us calculate an explicit expression for Cn,k- First, by the 
definition of A„.fc and A„ it is observed that c„.„ — 1, i.e. 
the orthogonal polynomials ^ are monicQ Furthermore, after 
inserting (|2]i into A„.fc and denoting q = e"^ , we obtain (see 
Appendix for details) 



A. 



E{n] 



^»nr=o^n"^»+i(g^-g') 
^^n;=.+i('z^"-'?'=)n-=o('?'=-9') 



where E{n) = e»("-i)M+V"(2»'-3n+i). por k = n, A^ ^ 
An n and ^ becomes 



ra-2 n-1 

A„=£;(n)n n (^^ 



(9) 



^Monic polynomial is defined as the polynomial where the coefficient of 
highest degree term is unity. 



After substituting (O and ^ into © we find that 



= ('—1 V"+fec("-'^')A'„("-l/2)(n-fc) 



c„,fc = (-l)'''^'-'e-'^ ---g^ 



(10) 



- — f. '- L,! — c^rri — ?T — ^ — - is the generalized 
binomial coefficient. 



where E]^ 



III. Approximation to the Density of Product of 
Nakagami-to Random Variables 

A. General Framework 

According to the CLT, the distribution of a product of 
RVs can be approximated by a lognormal distribution when 
the number of RVs is large. Motivated by this, we choose 
/ln(2:) as an initial approximation for the distribution of a 
product of RVs in the context of the moment based density 
approximation, which is derived in ifTOl . Let M{k) denote the 
fc-th moment of a density function f{x), the approximation to 
f{x) reads 



AT 



I{x) ^ /LN(a;) 2J m'^t{x) 

4=0 



(11) 



where the constant rji — jj- X]I=o Ci.kM[k). Note that the first 
TV moments of the approximation (fTTT i are matched with the 
corresponding moments of f{x) [10|. 

Alternatively, (fTTT i can be rearranged by combining the 
coefficients of 7ri(a;) with the same power of x, resulting in 
f{x) ~ hAx)Y.^^^£,,x' where ^^ = J2k=j^kjVk- After 
direct integration the approximated CDF attains the form 



N 



Fix)^J2^^''^^ 



log(a:) - fi 



(12) 



where $(•) is the CDF of a standard normal RV The CCDF 
F{x) = 1 — F{x) is approximated by replacing F{x) with the 
expression (fT2l i. 

It is a typical situation that the moments of certain RV 
are relatively easy to obtain whereas its exact distribution 
is difficult to calculate or unavailable [13|. The proposed 
approximations (fTTT i and (fT2l i are particularly useful in this 
situation, as will be shown in the next subsection. 

B. Moment Based Approximation to Distribution of Product 
of Nakagami-m RVs 

Let RV P — Y[i=i ^i ^^ ^ product of K Nakagami-yn RVs 
each with the PDF 

„2mi — 1 



^8) JbA^)-- 



2mT'x' 



nTT{m,) 



-exp 



^r 



x£ (0,00). (13) 



We first consider the case when the variables Ri are corre- 
lated with each other In the literature, there exists different 
representations for the joint PDF of multivariate Nakagami-m 
RVs, which are either limited in parameter values or cross- 
correlation structures. In this paper, we adopt a recent result 
derived in 1 14], which gives the joint PDF as a single integral. 
In the proposed approach, the RVs Ri are assumed to have 
the same fading parameter m and the joint PDF is valid 



for integer and half-integer values of m. The power cross- 
correlation coefficient between Rif and R^ is of the form 



E[flffl|]-E[i?j]E[J?|] 

PRf,R'i = ; = — = \^r 

Var[i?2]Var[i?2] 



(14) 



Note that the derived approximation framework is not limited 
by the above multivariate model. 

Based on the joint PDF 1 14, eq. (20)], the fc-th moment of 
the RV P is calculated by using [IS] eq. (6.643/2)] as 



K 



Mik)^ 3z.t'!lL um^~m^ 



TO^'=/2r(m)^+i 



K 



f^-^e-^ 



RiF. 



XH 



1 -t;^"^^ 



\f 



1 



dt, (15) 



where iFi(a, 6, c) denotes the Kummer confluent hypergeo- 
metric function. The parameters /i and a^ of /ln(2;) can be 
obtained by equating them with the mean and variance of 
log(P) respectively as 



M = f^E[log(i?0] = if^ 



1=1 

K 



i=l 



^oim) -log 



n, 



(16) 



a2=^ Var[log(i?,)] + 2^Cov[log(i?0,log(i?,)] 

i<j 
i<j 



i=l 
K 



i=l i<j ^-^0 



4r(m)e* 



d^-OO ,(17) 



where ^„(a;) is the polygamma function that is defined as 
the (n + l)-th derivative of the logarithm of r(x), and Q 
i ' 

2 



\l/o(™) - log ( rr) ■ ^^^^ *^ function Ii{t) is 
U,(l-Af) 



/,(t) = *oM+log 



iFi a, TO 



xn 



Xi-1 



where iFi (a, b, c) refers to the derivative of iFi(a, b, c) with 
respect to the parameter a. 

Next consider the case where the RVs i?, are independent. 
We notice that with A^ = 0, the joint PDF IHl eq. (20)] 
is in the form of a product of individual PDFs of Ri (i ~ 
1, 2, • • • , K), which implies statistical independence. Thus, we 
let Ai = and calculate the moments of P from ( fTSl l 



«(^)^n2^'^"-^" 



(18) 



While the calculation of /i is not affected by the independency 
and it is given by expression ( fT6] l. the variance a^ is reduced 
to the form 

1 ^ 
a'-^^^iM (19) 

1=1 

with all the covariances Cov[log(_Ri),log(i?j)] in ([IT) equal 
to zero. It is noted that ( fTSl l and ( fT9] l agree with the results 
given by [ 1] . 




Fig. 1. CCDF of products of six Nakagami-rra RVs when f2 = 1 and 
N = 16. (a) m = 1; (b) m = 4. Solid lines: approximative PDFs; markers: 
simulated PDFs. 



IV. Numerical Results 

In Fig. [T] we compare simulations with the approximative 
CCDFs for the product of six Nakagami-m RVs with pa- 
rameters m — 1, 4 and = 1, using orthogonal polyno- 
mials of up to 16-th degree. Here we consider equal cross- 
correlations, i.e. pjf! ^2 = p (i ^ j), and calculate the 
parameters of the approximation using (ITSt . (fT6t and (ITTt 
with p = 0.1, 0.5 and 0.8 for correlated Nakagami-m RVs. 
In case of independent RVs (with p = 0), equations ( fTSl l. (fTSl l 
and (fT9] l are applied. When m ~ 1 the RVs Ri are Rayleigh 
distributed, while with m ~ A the Ri are approximately Rician 
distributed with the Rician K-factor given hy k, — 6.46 [5 1. For 
each simulated CCDF curve, we generate 10^ realizations of 
the Nakagami-m RVs i?i (i = 1, • • • ,6) using Sim's method 
lfT6i . As a comparison, we also plot the approximative CCDFs 
calculated from [8] eq. (36)] for the cases with independent 
Nakagami-m RVs. Fig. [TJb) shows that differences between 
the approximative CCDFs and simulations are less than 10^'^ 
when /O = and 0.1, and less than 10"^ when p ~ 0.5. In 
addition, when m = 4 and p = 0, the proposed approximation 
yields improved accuracy compared with the results given 
in Is). However, as m — 1, the approximations with the 
same correlation coefficients are noticeably inaccurate due to 
slow convergence of the series in (fT2] i. Note that in both 
cases the approximative CCDFs with p = 0.8 deviate from 
the simulations since the CLT fails under high correlation 
condition. 

Table |I] summarizes the MSB e^ ^ ^^\F*{x) - 
F{x)\'^ dF{x) between the approximation F{x) and the empir- 
ical distribution function F* [x) of the product P as a function 
of the number of RVs K. When K is smaU (K < 10), 
Table U shows that the approximated CDF achieves consider- 
ably smaller MSB in low-correlation cases [p — 0.1 and 0.5) 
compared with the independent case (p — 0). As K increases, 
the accuracy for the independent case becomes superior to 
those of correlated ones. When m — 4, all approximations 
are much more accurate compared with the corresponding 
Rayleigh cases (m = 1), and the MSBs e^, except for those 



TABLE I 
Mean square error <? of approximative distribution 



K 


2 


4 


6 


8 


10 


12 


14 


16 


18 


20 


m = 1, p = 


1.14e-3 


1.09e-3 


6.28e-4 


3.78e-4 


2.72e-4 


2.02e-4 


1.67e-4 


1.45e-4 


1.28e-4 


1.14e-4 


m = 1, p = 0.1 


1.06e-3 


5.00e-4 


1.05e-4 


1.66e-5 


1.73e-5 


5.35e-5 


1.03e-4 


1.72e-4 


2.40e-4 


3.21e-4 


m = 1, p = 0.5 


9.99e-4 


1.68e-4 


2.10e-5 


9.54e-5 


2.07e-4 


3.18e-4 


4.16e-4 


5.36e-4 


6.24e-4 


6.98e-4 


m = 1, p = 0.8 


2.17e-3 


9.16e-4 


7.72e-4 


7.41e-4 


7.28e-4 


7.21e-4 


7.09e-4 


7.15e-4 


7.13e-4 


7.16e-4 


m = 4, p = 


8.13e-6 


2.29e-5 


2.43e-5 


1.64e-5 


2.24e-5 


3.15e-5 


4.34e-5 


4.78e-5 


5.14e-5 


5.02e-5 


m = 4, p = 0.1 


8.31e-6 


1.16e-5 


2.65e-6 


2.33e-6 


1.77e-5 


4.95e-5 


8.93e-5 


1.15e-4 


1.31e-4 


1.47e-4 


m = 4, p = 0.5 


7.03e-6 


1.34e-5 


7.44e-6 


9.73e-6 


1.74e-5 


2.48e-5 


3.36e-5 


4.21e-5 


4.99e-5 


5.49e-5 


m = 4, p = 0.8 


3.58e-5 


3.20e-4 


4.63e-4 


3.03e-4 


2.27e-4 


2.12e-4 


2.09e-4 


2.11e-4 


2.09e-4 


2.06e-4 



where X > 16 and p — 0.1, is less than 10^* in case of 
independent RVs or when correlation is low. 

Although the accuracies of the approximations (flU and (fTZt 
should be improved as N increases, analysis of the exact 
improvement is difficult due to the multiple nested finite 
summations involved. However, numerical results indicate that 
the «-th term in the summation (fT2] l converges quickly to 
zero since the coefficient ^iVi converges to zero and the CDF 
$ ( °st^^~M _ j(j J is bounded. The value A^ = 16 is shown 
to be sufficient to yield a stable CDF approximation. 

V. Conclusion 

Knowledge of the distributions of product of random vari- 
ables is important for understanding the performance of vari- 
ous communication systems. In this work, we first derived a 
closed-form expression for the orthogonal polynomials asso- 
ciated with general lognormal density. The derived result was 
subsequently applied in approximating the distribution of the 
product of random variables. As an example, we calculated 
closed-form approximations for the distributions of product 
of Nakagami-rn variates. Under light fading conditions, the 
resulting expressions achieve a good trade-off between com- 
putation complexity and approximation accuracy for small 
cross-correlations. It is still an open problem to find a closed- 
form (approximated) distribution for the product of correlated 
Nakagami-m RVs under deep fading condition. 

Appendix 
We note that the determinant 



,(»+j)m+- 



(20) 



i— O,--- ,n— 1; J— O,--- ,n; j^k 

the term e'*^"^' 

2„2 , 



remains unchanged by factoring out the term e*'^"'"* '^ /^ from 
the i-th row and the term e^^+-'''^ 1"^ from the j-th column. 
Namely, we have 



Vk 



i—O.--- ,n— 1; j—0,--- ,n; j^k 
n 



^nllLo lVi=^+l(<l^ ^ 1") 



'^kU'jLk+li'l' ^ 1'')Il^=oi<l 



,fe 



(21) 
(22) 
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equality follows from the fact that 
determinant. 
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